Hypothesis of strong chaos and anomalous diffusion in coupled symplectic maps 
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We investigate the high-dimensional Hamiltonian chaotic dynamics in A'^ coupled area-preserving 
maps. We show the existence of an enhanced trapping regime caused by trajectories perform- 
ing a random walk inside the area corresponding to regular islands of the uncoupled maps. As 
a consequence, we observe long intermediate regimes of power-law decay of the recurrence time 
statistics (with exponent 7 = 0.5) and of ballistic motion. The asymptotic decay of correlations 
and anomalous diffusion depend on the stickiness of the A'^-dimensional invariant tori. Detailed 
numerical simulations show weaker stickiness for increasing A'^ suggesting that such paradigmatic 
class of Hamiltonian systems asymptotically fulfill the demands of the usual hypotheses of strong 
chaos. 
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I. INTRODUCTION 

Hypotheses of strong chaos are keystone requisites for 
both basic theories (statistical mechanics ELQ) and com- 
putation methods (transport properties elimination 
of chaotic variables In all cases, the hypotheses of 
strong chaos consist essentially of two main assumptions 
about the (high-dimensional Hamiltonian) dynamics: (i) 
ergodicity, i.e., the existence of a single chaotic com- 
ponent and negligible measure of the regions of quasi- 
periodic motion; and (ii) strong mixing (exponential de- 
cay of correlations). Generic low-dimensional Hamilto- 
nian systems violate both hypotheses since the chaotic 
trajectories stick to the border of islands of regular mo- 
tion introducing long-term correlations in the dynam- 
ics 

For many higher-dimensional systems it is generally 
expected that hypothesis (i) is effectively valid since the 
measure of the quasi-periodic trajectories typically de- 
creases (exponentially) with the number N of degrees 
of freedom (see, e.g.. Sec. 6.5 of ref. [3])- However, 
the vanishing measure of regular regions [hypothesis (i)] 
does not guarantee the fast decay of correlations [hy- 
pothesis (ii)]. The breakdown of hypothesis (ii) occurs 
due to the trapping (stickiness) of chaotic trajectories 
around non-hyperbolic structures in the phase-space, and 
even zero measure sets (e.g., the bouncing ball orbits in 
the Bunimovich stadium billiard) can be responsible for 
the anomalous decay of correlations [1, Q . Considerable 
progress has been achieved for the problem of trapping 
in area- preserving maps {N — \) B, [Ei B 1 while only few 
numerical results are known for iV = 2-5 UilhJ. On the 
other hand, higher-dimensional systems are usually con- 
sidered for the calculation of different properties, such as 
relaxation phenomena ^12'] and diffusion [13, 15]. 

In this Letter we perform a detailed study of the trap- 
ping properties of Hamiltonian systems with an increas- 
ing number of degrees of freedom by coupling N area- 
preserving maps [13, [13] ■ This system typically has in- 



variant tori, but hypothesis (i) is satisfied for large N in 
the sense mentioned above (our results do not apply when 
regular behavior prevails for N ^ 00 [T^). For weak 
coupling we observe an intermediate regime of enhanced 
trapping which we show to exist also for area-preserving 
maps perturbed by noise. Asymptotically, our numerical 
results show that the trapping decreases with N indi- 
cating the effective validity of hypothesis (ii) for high- 
dimensional Hamiltonian systems composed of coupled 
low-dimensional ones. Additionally, we show how the 
different trapping regimes impact on the anomalous dif- 
fusion of the perturbed standard map and lead to a non- 
trivial dependence of the asymptotic diffusion coefficient 
on the perturbation st reng th, clarifying previous conflict- 
ing results [H [H [M [ia ■ 



II. COUPLED SYMPLECTIC MAPS 

We construct a time-discrete 2A^-dimensional Hamil- 
tonian system by the composition T o M of the indepen- 
dent one-step iteration of N symplectic 2-dimensional 
maps M = [Ml, . . . ,Mm) and a symplectic coupling 
T — (Ti, . . . ,T/v). As a representative example of 2-d 
maps we choose for our numerical investigation the stan- 
dard map : 



Pi + KiSa\{2'Kqi) mod 1 

q_i + Pi + Ki sin(27rgi) mod 1 



(1) 



and a coupling potential between the maps i,j given 
by Vij = cos[27r((7j — (7i)][22]. The action of the cou- 
pling on the i-th map is hence given by 



T,, 



E^i^ijSin[27r(gi-gj)] 



(2) 



which corresponds to a perturbation Api. The full cou- 
pling T is symplectic provided Ci.j = For simplicity 
we use all-to-all coupling with i^^.j = — in which 
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case a numerically convenient mean-field representation 
can be written ;12| 

^Pt = ^ sin[27r((7, - q/)] = ^\m\ sin(27r(?j - cj)), 

(3) 

where m = {m^,my) = ^y^J2j[cos{qj),am{qj)] 

and tan((/)) = my/mx- When the isolated systems are 
chaotic and weakly coupled we can assume that the po- 
sitions qj are uncorrelated and approximate each term of 
the sum in Eq. ([3]) by a random variable y e] — 1, 1[ dis- 
tributed according to the density = l/(27r-\/l — y"^)^ 
which has variance Oy = a/2/2. Additionally, for large N 
Eq. ([3]) tends to a normal distribution with a = £,ay, i.e., 
a finite perturbation strength. This is the main motiva- 
tion for the choice of the rescale factor I/^JN — 1 what is 
a major difference from other coupled-oscillators models, 
where a factor 1 /N is used instead [l^, [31 • 

III. TRAPPING IN THE MODEL 

Motivated by the previous considerations we study 
initially the trapping properties of one standard map 
{N = 1) perturbed by white noise, i.e, we replace Eq. ([3]) 
by Api = ^(5, where (5 is a Gaussian distributed ran- 
dom number with zero mean and variance a = ay. We 
fix Ki = 0.52 in ([1]), where a single large regular island 
is visible in the phase space (inset of fig. [1]), what fa- 
cilitates the interpretation of our results. The trapping 
is measured by regimes of power-law decay of the recur- 
rence time statistics (RTS) p(r) ~ , defined as the 
probability of a trajectory to return at a time T > r to 
a pre-defined region. During the long recurrence times 
the trajectory performs almost a quasi-periodic motion 
leading to a power-law decay of correlations with expo- 
nent 7c = 7 — 1 Hi . We have iterated a single trajec- 
tory 10^^ times and recorded the times T between suc- 
cessive recurrences to a large region away from the is- 
land [23]. The RTS for different noise intensities ^ is 
shown in fig. [TJ Three different regimes can be identified 
for small ^: 

(Rl) For short times (r < ti,2) the RTS follows the 
unperturbed one = 0), i.e., it shows an exponential 
followed by a power-law decay (with exponent 7^1). 

(R2) For intermediate times (ri^2 < t < T2^z) the RTS 
shows an enhanced trapping due to trajectories that en- 
tered the island through the action of the noise. Once 
inside the island the trajectories revolve around the cen- 
tral elliptic periodic orbit and perform a random walk 
in the perpendicular direction. The power-law exponent 
tends to the value of a random walker 7^2 ~ IRW = 0.5. 
As in typical RTS of Hamiltonian systems, we observe 
additional oscillations. 

(R3) For long times (t > T2,3) the RTS decays expo- 
nentially. 

We obtain now the dependence of ri,2 and T2,3 on ^. 
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FIG. 1: (Color online) (a) RTS for the standard map with 
K = 0.52. The thick solid line corresponds to the unper- 
turbed case (C = 0), whose phase space is shown in the in- 
set. The thin solid lines correspond to perturbations due to 
white noise with ^ = lO", 10~^ (from left to right). In 
each curve ri,2 (•) and r2,3 (■) are indicated. The dashed 
lines correspond to the cases of two (upper) and three (lower) 
coupled standard-maps with ^ = lQ~^,Ki — 0.5214, and 
K2 = 0.5108 pj. (b) Dependence of ri,2 and r2,3 on ^: 
ri,2 ~ 0.585"° (lower line) in agreement with Eq. (U); and 
T2,'i ~ 0.04^^^ (upper line) in agreement with Eq. ([5]). 

The starting time of (R2), ri,2, occurs roughly when the 
displacement due to the noise is of the same size of the 
distance to the regular island. Following the arguments 
of ref. one obtains 

ri,2 - (4) 

Numerical simulations in different maps show that /? ~ 1 . 
In ref. [lB| the noise intensity was compared to the size of 
the chaotic layer between two cantori close to the island. 
Using the Markov-tree model for stickiness introduced in 
ref. f5i| it was obtained that [3 — 1/(277^1 — 1). Since usu- 
ally 1 < jfji < 2 one obtains 1/3 < (3 < I. Differently 
from ref. 16], we do not observe an exponential decay 
right after this time. As argued above, the exponential 
regime (R3) appears at a later time T2^3 > T1.2 due to 
the finiteness of the random walk domain. Considering 
the measure /j-i of the (largest) island inside which a ran- 
dom walker (with step size proportional to 1^) performs 
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a diffusive motion, the dependence of T2,3 on ^ can be 
estimated as 



(5) 



In fig. [T]3 we verify tfie agreement of Eqs. ^ and ([5]) 
with the values estimated numerically. The numerical 
values of ti^2 and r2,3 were estimated as the intersecting 
point of two fitted power-laws (or exponential), which 
extend over more than two decades for small ^ and can be 
extrapolated for higher ^. The noise has two qualitatively 
different effects: while (R3) represents the typical cut- 
off of the power-law distribution [l^, during the novel 
regime (R2) the noise acts constructively (increasing the 
regularity of the dynamics) by allowing trajectories to 
penetrate regular islands. 

We investigate now the fully deterministic system 
given by the composition of Eqs. ([1]) and ([2]). The ex- 
istence of A'^-dimensional invariant tori is confirmed for 
small coupling by the application of the Kolmogorov- 
Arnold-Moser theorem to the tori built as a direct prod- 
uct of the 1-dimensional tori of the N uncoupled maps 
(the quasi-periodic orbits inside the island). All our 
numerical simulations for ^ < 0.1 are consistent with 
the assumption that these are the only existent tori, in 
agreement with a more general picture of Hamiltonian 
systems: exponential decay of the measure of the tori 
with TV {utori = l^'i'^ for ^ ^ 0) and nonexistence of tori 
of dimension smaller than N - "Froeschle conjecture" 0] . 
Furthermore, due to Arnold diffusion, it is consistent to 
assume that a single chaotic ergodic component exists. 
This means that the motion in the 2A^-dimensional phase 
space belongs to the chaotic component if the variables of 
at least one map belong to its chaotic component when 
uncoupled. 

The numerical RTS obtained for > 1 and for the 
noisy model are almost indistinguishable for recurrence 
times belonging to the regimes (Rl) and (R2) described 
above. Two representative examples - for ^ = 10"'^ 
and A^ = 2, 3 - are depicted as dashed lines in fig. [T(23|. 
Remarkable differences are observed in regime (R3). The 
detailed analysis shown in figH^ for ^ = 0.05 indicates 
that: 

(R3') For large times (r > T2^3) the RTS shows an 
exponential followed by a power-law decay (with expo- 
nent 7i?3) due to the stickiness to A^-dimensional tori. 

We investigate next the dependence on A^ of the 
asymptotic power-law exponent 7^3, which corresponds 
to a pure high-dimensional effect. The term asymptotic 
here has to be taken with caution since already for A = 1 
the convergence of the RTS to a well defined power- law is 
very slow in r [J]. However, for our purposes it is enough 
to perform a comparative analysis for different A^ and 
similar recurrence times and initial conditions. It is nec- 
essary to distinguish between two effects of A on the RTS 
seen in fig. [2^: the later onset of the power-law regime for 
increasing A^, related to the smaller measure of the tori ~ 
hypothesis (i); and the different values of 7/^3 (slopes of 
the tails), related to the stickiness to higher-dimensional 
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FIG. 2; (a) RTS for ^ = 0.05 and A = 2, 3, 4, 5, and noise 
(from top to bottom), (b) Survival probability f{t) of 10^*^ 
trajectories near the A-dimensional tori for ^ = 0.05 and 
N = 2, 15 (from top to bottom), (c) Exponent 7^3 fitted to 
the power-law regime of f{t) for ^ — 0.1, 0.05, 0.03 (from top 
to bottom). In all cases Ki = 0.52 and Ki G [0.51, 0.53], « G 
{2, N} were used in map ([l}. 



tori ~ hypothesis (ii) investigated here. Figurej2k already 
suggests that 77^3 increases with A'. In ref. similar 
results were reported for a different system and N = 2 
and 3. For improved statistics in the tails we study the 
survival probability f{t) inside a region containing the 
A-dimensional tori of 10^° trajectories started close to 
them 25]. f{t) is shown in fig.l^b where it is evident that 
the power-law exponent increases with A^. The estimated 
exponents shown in fig. [2j3 for different moderate values 
of the coupling ^ suggest a linear dependence oc A^. 
This result is consistent with hypothesis (ii) for high- 
dimensional systems since the sharp tails indicate fast 
decay of correlations. 



IV. ANOMALOUS DIFFUSION 

One of the most important effects induced by the 
stickiness is the anomalous diffusion of chaotic trajec- 
tories. We consider next the influence of the trapping 
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regimes (Rl-3) discussed above on the widely studied 
case of diffusion in the momentum p of the standard 
map d, 0, Ell- The diffusion is measured considering 
the map ([J) opened in p and calculating the temporal 
evolution of the dispersion of an ensemble of trajecto- 
ries: (Ap^(<)) = Df^ . Diffusion is normal \i v = 1. It 
is well-known that anomalous superdiffusion 1 < v < 2 
is obtained for the parameters K of the standard map 
where the so-called accelerator modes (ballistic islands 
of the open map) exist [1, 0, [ll]- The diffusion expo- 
nent u is directly linked to the exponent 7 (of the RTS 
due to the trapping around accelerator modes) by [31 



v — 



- 7 



if 
if 
if 



7 < 1 , 
1 < 7 < 
7 > 2. 



(6) 



Previous publications report numerical results that 
emphasize different effects of the noise on the diffusion: 
while in refs. |15|.ll6l] the onset of normal diffusion was ob- 
tained, in ref. [17| the possible enhancement of anomalous 
diffusion was reported. Applying the results described 
above we show that actually both effects exist for differ- 
ent time scales. 

The three trapping regimes (Rl-3) induce different 
diffusion regimes, as shown in fig. [3] and described be- 
low. During (Rl) the anomalous diffusion is similar to 
the unperturbed case. In (R2) the trapping exponent 
tends to 7^2 = 0.5 and ballistic motion = 2 is pre- 
dicted according to Eq. We have verified numer- 
ically that the beginning and end of this regime oc- 
cur at times proportional to, but greater than, those 
of (R2) (we denote these times by 2 and 3, respec- 
tively). Asymptotically the exponential RTS in (R3) im- 
plies normal diffusion [2^. The asymptotic diffusion co- 
efficient Da = limf^oo(Ap2(i))/i can be determined by 
the intermediate anomalous regimes and shows a nontriv- 
ial dependence on the noise intensity ^: for weak noise, 
Eqs. (dl and ([5]) indicate that the dominant contribution 
comes from the ballistic regime associated to (R2) and 

thus Da ~ [DB^iTlz — ^1.2)] ^ ^ For stronger noise 
the major contribution is given by a regime of superdif- 
fusion corresponding approximately to (Rl), that can be 
estimated as Da ~ {Drit\ 2)"^^ ^ ^-0(1^-1)^ where v is 
the unperturbed anomalous diffusion exponent (for small 
times). Considering the composition of this two effects 
and joining the multiplicative terms in two fitting param- 
eters a, h we obtain 



DA{0=a^ 



-Si{v-\) 



(7) 



In the inset of fig. [3] we show the remarkable agreement 
of the numerically obtained diffusion coefficient and ex- 
pression ([7]). Specially interesting is the transition for 
small ^ to an asymptotic ^"^ dependence of Da, which 
can be seen as a direct consequence of the nontrivial 
trapping regime (R2), is absent in the case of 1-d maps 
(e.g., Pomeau-Manneville maps) [l^ [l^, and was not 



previously reported in refs. [l5|, [iTj]. The asymp- 
totic scale was predicted through different arguments 
in ref. [l^l , and can be traced back to the works of Tay- 
lor M. 
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FIG. 3: (Color online) Diffusion in a standard map with K\ = 
1.07 (accelerator mode) coupled to noise (thin solid lines) and 
to another standard map with K2 = 0.52 (dashed lines). 
The uncoupled case is depicted as a thick solid line. From 
top to bottom i = 0.0003,0.001,0.005,0.1. 10** trajectories 
were used with initial conditions away from islands. When 
the curves become constant (for t — > 00) diffusion is nor- 
mal. Inset: asymptotic diffusion coefficient for different 
noise intensities. The solid line corresponds to Eq. ((T} with 
a = 0.604, h = 1.65 • 10"^ and -/3(1 - v) = 0.273. 

The diffusion theory described above is valid for noise 
perturbation and for deterministic high-dimensional sys- 
tems whenever 7^3 > 2 [see Eq. dS])]. When ^r-^ < 2 one 
observes additionally an asymptotic regime of anoma- 
lous diffusion, as shown in fig. ^ for N ~ 2 coupled 
standard maps and different coupling strengths ^. We 
see that the question of whether the asymptotic diffu- 
sion is anomalous or normal depends crucially on how 
the chaotic trajectories stick to A^-dimensional tori (ex- 
ponent 7i?3). Our results indicate that in general jr3 
increases with N (fig.J^l and the diffusion is normal, in 
agreement with refs. [isl. [T5|. However, our explanation 
for this behavior is not the absence of hierarchical tori or 
the Arnold diffusion, as argued in ref. [l^ for a similar 
system, but simply that 7/^3 > 2 for sufficiently large N. 



V. CONCLUSION 

We have shown that correlations decay faster as the 
dimensionality iV of a paradigmatic class of Hamiltonian 
systems increases. More precisely, our numerical sim- 
ulations show that the asymptotic power-law exponent 
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of the RTS increases with N, providing evidence that, 
for large enough N and times, coupled symplectic maps 
can be considered for all practical purposes as ergodic 
and strongly chaotic. This suggests a novel explanation 
for the onset of strong chaos in high-dimensional Hamil- 
tonian systems whose generality has to be determined 
through the specific study of other classes of systems. An 
unexpected behavior is that for small coupling strength 
or noise intensity a long intermediate enhanced-trapping 
regime with 7^2 ~ 0.5 exists due to the trapping in- 
side remnants of lower-dimensional regular regions. This 
implies a regime of enhanced anomalous diffusion and 
a non-trivial dependence of the asymptotic diffusion co- 
efficient on the perturbation strength, given by Eq. ([7])- 
Applications of these results include noise-perturbed low- 



dimensional system as well as fully deterministic ex- 
amples such as galaxy dynamics [ll| and active trans- 
port (I'sj . Further investigations are needed to determine 
whether our results provide a valid description also of sys- 
tems like Hamiltonian mean field models, where similar 
transient regimes were observed p^ . 
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